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ABSTRACT 

Spectroscopic analyses of gravity-sensitive line strengths give growing evidence towards 
an excess of low-mass stars in massive early-type galaxies (ETGs). Such a scenario 
requires a bottom-heavy initial mass function (IMF). However, strong constraints can 
be imposed if we take into account galactic chemical enrichment. We extend the anal¬ 
ysis of Weidner et al. and consider the functional form of bottom-heavy IMFs used 
in recent works, where the high-mass end slope is kept fixed to the Salpeter value, 
and a free parameter is introduced to describe the slope at stellar masses below some 
pivot mass scale (M<Mp = 0.5Mq). We find that no such time-independent parame- 
terisation is capable to reproduce the full set of constraints in the stellar populations 
of massive FTGs - resting on the assumption that the analysis of gravity-sensitive 
line strengths leads to a mass fraction at birth in stars with mass M < Q.5Mq above 
60%. Most notably, the large amount of metal-poor gas locked in low-mass stars dur¬ 
ing the early, strong phases of star formation results in average stellar metallicities 
[M/H]< —0.6, well below the solar value. The conclusions are unchanged if either the 
low-mass end cutoff, or the pivot mass are left as free parameters, strengthening the 
case for a time-dependent IMF. 
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1 INTRODUCTION 


The distribution of stellar masses in galaxies, i.e. the stellar 
initial mass function (IMF), is of fundamental importance 
for describing stellar populations. Amongst other things, the 
IMF drives the chemical evolution of galaxies and defines 
the stellar mass-to-light ratios (M*/L). The IMF has been 
intensely studied in the Milky Way (MW) and nearby dwarf 
galaxies, and it has been foun d mostly invariant for a large 
range of physical p arameters (iKroupa l2002l : IChabrieill2003l : 
iKroupa et M.l[201^ . 

However, in recent years this universality of the 
IMF has been called into question by a number of 


new ODservaxionai results \see, e.g. 
Hoversten & GlazebrookI 20081: 

iL^enarro et ai.i 

Meurer et al.l 

zuuo:: 

200c; 

Gunawardhana et al.l 201ll: van Dokkum & Conrov 

20101. I2012I: iGaDoellari et al.l I2OI2I: 

Ferreras et al.l 20131: 


Weidner et al. 


20iA|. 


On a first glance, not all of these 
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results agree well w it h eac h other. While, for example, 
iGunawardhana et al.l ll201lli find that the IMF becomes 
increasingly top-heavy in galaxies with a high star for¬ 
mation rate, other studies suggest bottom-heavy IMFs 
for massive early galaxies (ETGs), systems which are 
believed to have formed in a massive starburst. More hints 
about top-heavy IMFs have also been foun d in the bulges 
of M31 and the MW llBallero et al.l l2007l i as well as in 
massive globular clusters and ultracompact d warf galaxies 
llDabrinAausen et al.l[201^ : iMarks et al.ll201^ ') . 

The rise of large and deep galaxy surveys in recent years 
opened a new angle on ste l lar po pulations in galaxies. For 
example, iGappellari et ah] (l2012l i used integral field spec¬ 
troscopy and photometry of a volume-limited catalogue of 
260 ETGs to constrain the populations in these galaxies. 
They found that, independent of the choice of dark mat¬ 
ter halo model for the galaxies, the SDSS r-band M/L ra¬ 
tios suggest a syste matic transition from a standard IMF 
(e.g. iKroupal 1200 il l towards a distribution with heavier 
stellar M/L in the more massive galaxies (by ~60% at 
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a ~ 300 kms“^). Bottom-heavy IMFs for ETGs had already 
been s uRgested by studies of abundance line indicators in 
ET Gs JC enarro et al.ll2 003h an d bulges of late-type galax¬ 
ies (iFalcon-Barroso et al.l2003l . The need for non-standard 
IMFs had also b een shown from abundanc e ratios and the 
colours of ETGs JVazdekis et a 1ll996l. ll997l i. More recently, 
Ivan Dokkum fc GonrovI (l2010l . 2012l l analysed the gravity- 
sensitive indices NaSlOO and FeHO.99 and found strong evi¬ 
dence for a bottom-heavy IMF in ETGs. 

This view is however not without challenge. 
ISmith fc LucevI (|2013l ~l showed a strong lens massive 
ETG and compared the stellar lensing mass with mass 
estimates from population synthesis modelling, They found 
no striking deviation from a standard IMF. This galaxy 
is, however , very extended wh ereas compactness has been 
suggested (|Lasker et al.l l2013h as a driver for the IMF 
changes. 


Galactic chemical enrichment provides an additional 
constraint on the IMF, as the distribution of stellar masses 
plays an essential role i n the e nrichment of stellar popula¬ 
tions. In IWeidner et al.l (l2013l 'l it was shown that a time- 
independent standard bimodal IMF -- with a power-law at 
the high ma ss end, and a smooth tapering at low masses, 
as defined in IVazdekis et al.l (|l99fih - was not capable of ex¬ 
plaining all the observables of massive ETGs. The assump¬ 
tion of a bottom-heavy IMF results in low overall metallic- 
ities for the stellar populations, in contrad iction with the 
observations . The solution, first proposed in Vazdekis et ^ 
lll996l. Il99'it) . and revisited in IWeidner et al.l ll2013ll in¬ 
volves a time-dependent IMF, where a top-heavy phase, ex¬ 
pected during the first stages of evolution, is followed by a 
bottom-heavy phase. This scenario is physically motivated 
by the fact that the strong starbursts expected in mas¬ 
sive galaxies at high redshift inject vast amounts of energy 
into the ISM, changing the physical conditions to a highly 
turbulent medium at high pressure, perhaps inducing a 
stron g fragmentation process (lHopkinj2013l : lGhabrier et al.l 
201^). Howeve r , ther e is an additional aspect not covered in 
Weidner et al.l (|2013h : the chosen bimodal IMF tightly links. 


by construction, the high-mass end to the low-mass end, so 
one could still envision a distributi on where the high-mass 
end is kept at a Salpeter-like value (ISalpeteilll955li . whereas 
the slope at the low-mass end is left as a free parameter. 
Such a functional form of the IMF is adopted by some of th e 
groups in this field (see, e.g. iGonrov fc van Dokkurnll2012h . 
In this letter, we explore the consequences of this approach 
from the point of view of galactic chemical enrichment, and 
we find a significant mismatch to explain simultaneously the 
age-, metallicity- and gravity-sensitive features of massive 
galaxies. Therefore, the need for a time-dependent IMF is 
more compelling. 


This letter is structured as follows: ^provides a generic 
working definition for the functional form of the IMF. In ij3] 
the chemical evolution model with the variable IMF used 
here is described. The results from our model calculations 
are presented and discussed in @ Finally, Sj5] summarizes 
the conclusions. For reference. Tab. [T] shows conservative 
constraints from the literature on the observables used to 
test the hypothesis of a time-independent IMF. 



Figure 1. Functional form of the Initial Mass Function adopted 
here, along with standard choices from the literature (see text, 
and equation 1 for details). 


2 A FUNCTIONAL FORM OF THE IMF 


We assume a simple description for the initial mass function, 
as a truncated power law, namel^Q (see Fig. [T]): 


^(m) = 


_ ^ _ J Nim 
\ Nsm 


dm 


-{H-ri) 

5 

-(i-l-ra) 


Mj^ow < m ^ Mp 
Mp <m< Mhigh, 


( 1 ) 

where the slopes are given with respect to log-mass units. 
We assume the high-mass end branch to have a Salpeter- 
like slope, i.e. F 3 = 1.3. Our fiducial case fixes the mass 
scale of the pivot point at Mp = 0.5 Mq; and the mass 
range: Mlow = O.IMq; Mhigh = 100 M(t^. N ote that the 
case Fi = 1.3 is similar to a Salpeterl lll955h IMF. Varia¬ 
tions of this expression have been used in the lite r ature to 
descr i be the IMF fsee, e.g., iMiller fc Scald Il979l : iKroupal 
I 2 OOII : iMarks et al.l I 2 OI 2 I ) and have been used to explain 
the sign ature of low-mass stars in t he spectra of ellipti¬ 
cals fe.g. lGonrov fc van Dokkuml20l3) as well as analyses of 
the u nderlying chemical enrichment associat ed to these sys¬ 
tems (IVazdekis et al.iri996l . ll9^ : lBekkl2013l ). Spectroscopic 
data from gravity-sensitive features can only constrain the 
mass fraction in low-m ass stars at birth (see Fig. 21 of 
iLa Barbera et al] l2013l ). leaving a significa nt uncertainty 
with respect to the stellar mass-to-light ratio (iFerreras et al.l 
I 2 OI 3 I) , and compl icating comparisons with dynam ically- 
based studies (e.g. ICappellari et al.l[201^ : ISmithli2014l ). The 
assumption of a functional form of the IMF is needed for an 

adequ ate comparison. _ 

In lLa Barbera et al.l ll2013h it was found that regardless 
of the choice of IMF, the gravity-sensitive spectral features 
imply a conservative constraint on the stellar mass fraction, 
at birth, in stars with M<0.5Mq (Fo.s). Massive ETGs are 


^ Not e we use the 1, 3 subindex notation as in, e.g. iKroupa et ^ 
ll2013h . The missing subindex 2 would refer to an intermediate 
mass region that we do not consider here. 
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Figure 2. Contours of Fq.s as a function of the two IMF slopes, 
Fi and Fs (see text for details). The thick vertical dashed line 
represents the fiducial set of models explored in this letter, along 
with a top- + bottom-heavy model (cross, see §3). The circle 
corresponds to a standard Salpeter-like IMF. Values Fq.s ^ 0.6 
are expected for massive ETGs (orange contours). 

found to have Fb.s ^ 0.6. With the assumptions adopted 
here, no model with Fi < 1.5 is capable of explaining the 
IMF-sensitive spectral features of massive ETGs (see Fig. 13). 


3 CHEMICAL ENRICHMENT MODELLING 

AND OBSERVATIONAL CONSTRAINTS 

We explore a set of phenomenological models tracking galac¬ 
tic chemical enrichment through a reduced set of parame¬ 
ters. T hese models are presented i n detail in iFerreras fc Silla 
ll2000al f9). In IWeidner et alJ (l2013l ~) we apply them to a time- 
dependent IMF in order to explain the gravity-sensitive line 
strengths found in massive ETGs. In a nutshell, the buildup 
of the stellar component of a galaxy is described by four 
parameters: a gas infall timescale (r/), a star formation ef¬ 
ficiency (Ceff), that follows a Schmidt law, a formation red- 
shift (zfor) at which the whole process starts, and a fraction 
of gas ejected in outflows (Bout). 

We ran a grid of models adju sted to the stellar popu - 
lations in massive ETGs (see, e.g. Ide La Rosa et al.ll201lli . 
In order to achieve a homogeneously old population, we 
need to assume an early start for the star formation process 
(^FOR = 3), and negligible outflows (Bout = 0). Ghanges 
in these two parameters will mostly induce an overall shift 
in the average age and metallicity, respectively. Further¬ 
more, non-neglible outflows would produce lower metallic- 
ities, hence strengthening our conclusion towards a time- 
dependent IMF (see §4). The other two parameters, namely 
the gas infall timescale (r/) and star formation efficiency 
(Geff) are left as free parameters in the grid. The model 
grids are run for a range of IMF slopes, 1 < Fi < 2.5. 

Fig. [3] shows the results of the chemical enrichment 
modelling for three choices of IMF slope: Salpeter (Fi = 
Fa = 1.3, left); bottom-heavy (Fi = 2.2, Fa = 1.3, mid¬ 
dle), and the additional case of a top -I- bottom heavy IMF 
(Fi = 3, Fa = 0.9, right). For the latter, we modify the 
high-mass slope, Fa, to values that would be compatible 
with the obs ervations of top-heavy IMF s in star-forming 
systems (e.g. iGunawardhana et aLdinil . changing in ad- 




Figure 3. Contour plots with the main observable constraints 
from the chemical enrichment code corresponding to a massive 
early-type galaxy (i.e. with a stellar mass > IO^^Mq). Three dif¬ 
ferent choices of IMF slope are considered, a Salpeter-like case 
(top), a bottom-heavy case (middle), and a top -|- bottom-heavy 
case (bottom). The contours that are compatible with the obser¬ 
vational constraints (Tab. (TJ are shown as thick orange lines, the 
remainder as dotted contours. The contours in age and tgF are 
spaced in intervals of 0.5 Gyr. The contour intervals for [M/H] and 
M(< Zq/10) are 0.05 dex and 5%, respectively. All the contours 
in M/L are labelled, and measured with respect to the SDSS-r 
band. See text, and Tab. [T] for details. 
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Observable 

Constraint 

Reference 

Age (Gyr) 

[8,10] 

(1) 

[M/H] 

[-0.1,-1-0.2] 

(1) 

tsF (Gyr) 

[0.5, 2.0] 

(2) 

Ms(< Zq/10) 

[0.05,0.20] 

(3) 

Fo.5 

[0.6, 0.8] 

(4) 

Tr/T,,Q 

< 7.0 

(4) 


Table 1. Constraints on the general properties of the unre¬ 
solved stellar populations in massive early-type galaxies. The 
uncertainties are rough esti mates, quoted at a c onservative Itr 
level. The r eferences are: ( l)lTra ger et al.l ||200(J): iThomas et al.l 
ll2005|b__f2j_[de_Ljrjtosaj1^1.l ll201lh : 131 IVazdekis et^i^ ~ i^l^7^ : 
f4l iLa Barbera*er'arni2013ll 

dition the low-mass slope, Fi, to accommodate the high 
values of F0.5. For each case, we show (counter-clockwise 
from the bottom-right) contours of stellar mass-to-light ra¬ 
tio (in the SDSS-r band); mass fraction in low-metallicity 
stars (M(< Zq/IO)); star formation lapse (tsp); average 
age, and average metallicity ([M/H]). Averages are mass- 
weighted. The star formation lapse is defined as the time 
period between the 25th and 75th percentiles of the stellar 
mass build-up in the galaxy. tsF is therefore tightly linked to 
[Mg/Fe]. As a very conservative estimate, we assume that 
values tsF -J 2 Gyr are in conflict with t he observed over¬ 
abun dance of [Mg/Fe] in massive ETGs llde La Rosa et ahl 
[20 i3). M(< Zq/IO) is defined as the fraction of stellar mass 
with metallicity below 1/10 of solar. This is a n indicator of 
the G-dwarf problem llPagel fc Patchet^ll975l i. It is a well- 
known fact that some models can produce old, and overall 
metal-rich populations, although with a significant tail of 
low-metallicity stars. Models with a mass fraction over 10- 
20% in low-metallicity stars wo uld be in conflict wi t h the 
observations of massive E TGs (|\5rzdekis_et_^ 1 19961 . 1 19971 : 
iMaraston fc Thomasl[200^ : iNolan et al.ll2007l i. 


4 RESULTS & DISCUSSION 

Fig. E] shows that the Salpeter model (left) is capable of 
recreating the old, metal rich, [Mg/Fe] over-abundant pop¬ 
ulations without a significant low-metallicity trail. However, 
the fraction in low-mass stars (Fb.s = 0.44) is in conflict 
with the recent interpretation o f gravity-sensitive spectral 
features dLa Barbera et m]| 2013I ~). The top- -I- bottom-heavy 
model (right) is challenged by the overly high values of 
M/L in the region of parameter space compatible with the 
age/metallicity constraints. The same result would hold if 
we chose to keep Fa as a free parameter. Therefore, we rule 
out the option of a change in both slopes, and hereafter fo¬ 
cus on the case where the high-mass end slope is fixed to 
the Salpeter-like value (F3 = 1.3). The bottom-heavy model 
(middle) allows for a higher low-mass fraction (F 0.5 = 0.65) 
at the price of locking too much gas in low-mass stars during 
the early (metal-poor) stages. The average metallicities are 
significantly lower than those derived from the observations 
in massive ETGs. 

One could consider additional changes of the functional 
form of the IMF, most notably changing either the pivot 
point (Mp) " motivated by a change in th e physical pro p- 
erties of the star forming regions (see, e.g. iLarso nl l2005li - 


or the low-mass end ( Mlqw) ~ as suggeste d in the recent 
analysis of ETG lenses (iBarnabe et ahluOPTll . Note that we 
emphasize in this letter that all these changes would relate 
to an otherwise time-independent IMF. However, from the 
previous figure, it is expected that such changes neverthe¬ 
less lock large masses of low-metallicity stars, leading to re¬ 
sults that are incompatible with the metal-rich populations 
found in massive galaxies. To further illustrate this point, 
we show in Fig [4] a estimator based on the observational 
constraints presented above (see Tab. [TJ, where Gaussian 
constraints are imposed at the 1 a level over the allowed 
intervals, in effect producing very conservative limits. For 
M/L - where the constraint derives from dynamical M/L 
measurements - we simply penalize the likelihood for values 
higher than Tr = 7Tr,0 using a Gaussian with a — 0.5. 
We consider several cases, where either a range for the low- 
mass end (Mlow, top panel) or the pivot mass (Mp, bottom 
panel) are explored in a time-independent IMF. The solid 
lines correspond to different cases where all four constraints 
are imposed. The dashed lines show, for comparison, the 
case where the F0.5 constraint is removed from the analysis. 
As reference, the horizontal dotted line in bot h panels repre¬ 
sent th e hducial, time-dependent models A of lWeidner et al.l 
ll2013ll . where a top-heavy IMF is followed by a sharp tran¬ 
sition towards a bottom-heavy IMF after 0.3 Gyr. Note the 
different behaviour with respect to changes of either the low- 
mass end of the IMF, or the pivot mass scale. In the top 
panel, a transition is apparent at Mlow = O.IOMq, above 
which very bottom-heavy IMFs seem to be favoured, with 
the best case at O.lhM©. Nevertheless, even this option is 
rejected with respect to a time-dependent scenario. 

Fig. [4] confirms that for a wide range of options, a 
time-independent IMF is incompatible with constraints from 
galactic c hemical enrichment. Although the proposal pre¬ 
sented in IWeidner et al.1 (l2013ll is a simple toy model, this 
letter supports the need for a time-dependent mechanism 
that tips an initial top-heavy IMF in a strongly star-forming 
system, towards a bottom-heavy IMF during the final stages 
of the starburst. Such a mechanism would reconcile the 
apparent contradiction between the observations of star¬ 
bursting systems and the properties of quiescent galaxies 
that underwent a star-bursting phase in the past. We em¬ 
phasize that the relatively long duration (~ 1 — 2 Gyr) of a 
strong star formation phase in massive ETGs at high red- 
shift is expected to cause this transition. 


5 CONCLUSIONS 

We explore a model of galactic chemical enrichment, with 
the assumption of a time-independent IMF, with several free 
parameters controlling the contribution of low- and high- 
mass stars. The functional form (presented in [J2} is repre¬ 
sentative of the typical functions explored in the literature. 
A comparison of our models with conservative observational 
constraints in massive ETGs (Tab.[T]) reject this hypothesis, 
mainly based on the locking of too many low-mass metal- 
poor stars during the Erst phases of formation, and on the 

^ We emphasize t hat th e functional form of the (bimodal) IMF 
in IWeidner et al.l 1 I 2 OI 3 II changes both low- and high-mass ends 
with a single parameter, F. 
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Figure 4. Fits to the extended models, where the low-mass end of 
the IMF (top), or the pivot mass (bottom) are considered as free 
parameters in a time-invariant IMF. The values are shown for 
several cases (solid lines). We illustrate the constraining power 
of Fq.s by including the case where this constraint is removed 
from the analysis (dashed lin es). For reference, th e x^ of ^ time- 
dependent IMF (model A in IWeidner et al1l2013l) is shown as a 
dotted horizontal line in both panels (where all constraints from 
Tab.^ including Fq.s, are imposed). 


assumption that the recent observations of gravity-sensitive 
line strengths result in a constraint on the mass fraction 
in low-mass stars at birth (F 0.5 > 0.6). The best-fit time- 
independent models give stellar metallicities [M/H]'^ —0.6, 
i.e. significantly lower than the observational constraints - 
although higher than the metallicites expec ted for a time- 
indep endent “bimodal” IMF as defined in IVazdekis et al.l 
(Il996t) . In contrast, a simple time-dependent model, justified 
by the large energy injection during the strong star bursting 
phase in massive galaxies, leads to a consistent picture with 
an overall old, metal -rich and bottom-heavy stel lar popula¬ 
tion, as sug gested by Vaz dekis et ^ (Il996l . 199li l and, more 
recently bv lweidner et ah 11201^). 
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